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We obtain a Bahadur representation for sample quantiles of nonlinear func- 
tional of Gaussian sequences with correlation function decreasing as fc~" for 
some Q > 0. This representation is derived under a mimimal assumption. 



1 Introduction 

We consider the problem of obtaining a Bahadur representation of sample quantiles in a 
certain dependence context. Before stating in what a Bahadur representation consists, 
let us specify some general notation. Given some random variable Y, F(-) = Fy{-) is 
referred as the cumulative distribution function of y, ^(p) = (p) for some < p < 1 as 
the quantile of order p. If F{-) is absolutely continuous with respect to Lebesgue measure, 
the probability density function is denoted by /(•) = /y(-). Based on the observation 
of a vector Y = (1^(1), . . . ,Y(n)) of n random variables distributed as Y, the sample 
cumulative distribution function and the sample quantile of order p are respectively 
denoted by Fy (•; Y) and ^y {p; Y) or simply by F{-; Y) and ^ (p; Y). 

Let Y = (Y{1), . . . , Y{n)) a vector of n i.i.d. random variables such that F"{(^{p)) 
exists and is bounded in a neighborhood of ^(p) and such that F'{^{p)) > 0, Bahadur 
proved that as n — > +00, 
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with r„ = Oa.s. log(n)'^/'^) where a sequence of random variables Un is said 

to be Oa.s. (vn) if Un/vn is almost surely bounded. Kiefer obtained the exact rate 
72~3/4 log log (n)^/^. Under an Assumption on F(-) which is quite similar to the one done 



mixing variables, and recently in ml (I2OO5I) for short- range and long- 
linear processes, following works of 



es ha ve been pur- 



by Bah adur, extensions of ab ove results to dependent ra ndom variab 
sued in 



Sen and GhoshI (jl972l ) for i nixing variables, in lYoshiharal (119951 ) for strongly 



r ange 



dependent 

Hessd (I1990I) and IHo and Hsing I (|l996l ^. Finally, 



such a representation has been obtained by lCoeurjollyi (|2007l ) for nonlinear functional of 
Ga ussian sequeii ces with correlation function decreasing as k~°' for some a > 0. 



Ghosh (1971 



proposed in the i.i.d. case a much simpler proof of Bahadur's result 
which suffices for many statistical applications. He established under a weaker assump- 
tion on F{-) (-F'(-) exists and is bounded in a neighborhood of ^(p) and f{^{p)) > 0) 
that the remainder term satisfies r„ = Op{n~^^'^), which means that n^/^r„ tends to 
in probability. This result is sufficient for example to establish a central limit theorem 
for the sample quantile. Our goal is to extend Ghosh's result to nonlinear functional of 
Gaussian sequences with correlation function decreasing as A;^". The Bahadur represen- 
tation is presented in Section [2] and is applied to a central limit theorem for the sample 
quantile. Proofs are deferred in Section [3l 



2 Main result 

Let {Y{i)}f^ be a stationary (centered) gaussian process with variance 1, and correla- 
tion function p{-) such that, as i — > +00 

\p(i)\ ~ i-'^ (1) 

for some a > 0. 

Let us recall some background on Hermite polynomials: the Hermite polynomials 
form an orthogonal system for the Gaussian measure and are in particular such that 
E (i/j(y)fffc(y)) = j! Jj fc, where Y is referred to a standard Gaussian variable. For 
some measurable function g{-) defined on M such that E(g(y)^) < +cxd, the following 
expansion holds 

9{t) = E 1 ^i*^ = E i9{Y)H,{Y)) , 
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where the integer r defined by r = inf {j > 0, Cj 7^ 0}, is called the Hermite rank of the 
function g. Note that this integer plays an important role. For example, it is related to 
the correlation of g{Yi) and g{Y2), for Yi and I2 two standard gaussian variables with 
correlation p, since 'E{g{Yi)g{Y2) = Y.k>T %r/''' = ^ (^^)- 

Our result is based on the assumption that i^^^y^(-) exists and is bounded in a neigh- 
borhood of_^XEl:_JIl3isj^_achieved_^^ g{-) satisfies the following assumption 
(see e.g. iDacunha-Castelle and Dufid (jl982l ). p. 33). 

Assumption A(^(p)) : there exist C/j, i = 1, . . . ,L, disjoint open sets such that 
Ui contains a unique solution to the equation g{t) = Cg(y)(p)) such that F'^^Y){i{p)) > 
and such that 5 is a C^— diffeomorphism on uf^iUi. 
Note that this assumption allows us to obtain 



is referred to the probability 



Fg{Y)iCg{Y){p)) = fg(Y){^giY){p)) = 

where gi{-) is the restriction of g{-) on Ui and where 
density function of a standard Gaussian variable. 

Now, define, for some real u, the function hu{-) by: 

Kit) = l{g{t)<«}(t) - Fg(^Y)iu)- (2) 

We denote by r(n) the Hermite rank of /lu(-). For the sake of simplicity, we set Tp = 
T{S,g(Y)ip))- For some function g{-) satisfying Assumption A(^(p)), we denote by 

Tp = inf r(7), (3) 

7eufLig{C/,) 

that is the minimal Hermite rank of hu{-) for u in a neighborhood of £,g(Y){p)- Denote 
also by Cj{u) the j-th Hermite coefficient of the function hu{-). 



Theorem 1 Under Assumption A(^(p)), the following result holds as n ^ +00 

P-F{^g^Y){p);9(Y)) 



Hp;9iY))-^g(Y)ip) 



+ Op (r„(a,rp)) 



(4) 



fgiY){CgiY){p)) 

where g{Y) = {g(Y{l), . . . ,g{Y{n))), fori = 1, . . . ,n and where the sequence {rn{a,Tp))^-^^ 
is defined by 

ifoTp > 1, 

ifaTp = l, (5) 
ifaTp < 1. 



rn{a,Tp) 



n 



1/2 

log(n)i/2 
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Remark 1 The sequence Vni^oi^Tp) is related to the behaviouv shoTt-vange ot long-range 
dependent behaviour of the sequence hu{Y{l)), . . . , hu{Y(n)) for u in a neighborhood of 
S,{p). More precisely, it corresponds to the asymptotic behaviour of the sequence 

1/2 



Corollary 2 Under Assumption A{^(p)), then the following convergence in distribu- 
tion hold as n ^ +oo 
(i) if OTp > 1 



(6) 



where 



a, 



" f{p) 



1 ^ C'(p)^ 

Yl Yl ~^P^^y ^^^^ f^P^ = fg{Y){ig{Y){p)) and Cj{p) = Cj{ig^Y){p))- 



i&L j>Tp 



(a) if aTp < 1 



({p;9iY))-CgiY){p) 



d c-^{p) 



Tplfip) 



(7) 



where 



Ztj, = K{Tp, a) 



i(Ai + ••• + AtJ 



and 



K{Tp,a) 



{1 - aTp/2){l - aTp) 
Tp! (2r(a) sin(7r(l - a)/2)y 



1/2 



The measure B is a Gaussian complex measure and the symbol J' means that the domain 
of integration excludes the hyperdiagonals {Aj = ±Aj, i ^ j}- 

The proof of this result is omitted since it is a direct application of Theorem [1] 
and general limit theorem s a dapted to nonlinear function al of Gaussian sequences, e.g. 



Breuer and Majoii (|l983l ) and iDehling and Taqqu I () 19891 ). 
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3 Proofs 



3.1 Auxiliary Lemma 

Lemma 3 For every j > 1 and for all positive sequence {un)n>i such that Un — ^ 0, as 
n — ^ +00, we have, under Assumption A(^(p)) 



Hj{t)(p{t)l{lg^t}-^^^Y){p})\<Un}dt ~ Un Kj, 



(8) 



where t\jj %s defined, for every j > l,by 



Kj = < 



if 3 = 1, 



2, lyy^L ^'^-^Kfr^(^(p)) 



(9) 



Proof. Under Assumption A(^(p)), there exists no G N such that for ah n > no, 

L 



I = Y^Ii with Ii= f Hj{t)cP{t)l 
1^1 JUi 



dt. 



(10) 



Assume without loss of generahty that the restriction of g{-) on Ui (denoted by gi{-)) is 
an increasing function, we have 



= / mm) 

JUi 

[■9i'^{iip)+u„) 

= / Hj{t)(p{t)dt 



h^{p)-un < g{t) < ^{p)+«„}^* 



4>{'mi,n) - 4>{Mi^n) = ("^^,n - Mi^n) if J = 1 

(</,0--i)(M,,„) - (t>^i-^\m,,n)) if i > 1, 

where Mj,„ = g~^{^(jp) + Un) and mi^n = g~^{^(jp) - Un). Then, there exists ujn,i,j e 
[mi^n, Mi^n] for every j > 1 such that 



Ii=< 



{rrii^n - Mi^n) <P^'^\^^n,i,l) if J = 1 

^ (-l)^' (Mi,„ - mi,n) 4'^'-''^ {^n,i,j) if i > 1. 

Under Assumption A(^(p)), we have, as n — +oo 

~ 9i^mP)) and Mi,„ - rui^n ~ 2«„ 

9'{9i (up))) 



which ends the proof. 
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3.2 Proof of Theorem [T] 

For the sake of simplicity, we set ^{p) = ^{p;g{Y)), ^(p) = Cg(y)(p), F (•) = F{-;g{Y)), 
H-) = Fg{Y){-) et /(•) = fg(Y){-) and r„ = r„(a,rp). Define, 

Vn = r-'(cip)-aP)) and = . 

P 

The result is established if 1^ — Wn — > as n +oo . It suffices to prove that Vn and 
Wn satisfy the conditions of Lemma 1 of iGhoshI (jl97ll ): 

• condition (a) : for all 6 > 0, there exists e = e{6) such that P {\Wn\ > e) < 6. 

• condition (b) : for all y G M and for all e > 

lim P (K <y,Wn>k + e) and lim P (K > y + e, VF„ > k) 

n— »+oo n— ++00 

condition (a) : from Bienayme-Tchebyshev's inequality it is sufficient to prove that 
EW2 = 0(1). Rewrite Wn = ^ EH^+i hip) (^(0)- Let cj (for some j > 0) denote 
the j-th Hermite coefficient of /i^(p)(-). Since ^5(p)(') has at least Hermite rank Tp, then 

E^" = ^ E E(/i^(,)(y(n))%rt(n^2))) 



il,«2=l 
2 " 



^ E E c,iC,,E(i/,,(y(n))^,,(ni2))) 

«1,*2 = 1 jl,32>Tp 



from Remark [TJ 

condition (b) : let y G M, we have 

{Vn < y} = {e (P) < 2/ X r„ + e(p)} 

= {p<F(2/xr„ + e(p))} ={^„<2/n}, (11) 

with 

= [F{y X rn + e(p)) -f(^+ ^p) 
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and 

,-1 



Under Assumption A(^(p)), we have ^ y, as ri +00. Now, prove that Zn—Wn 0. 
Without loss of generahty, assume y > 0. Then, we have 

-1 1 " 
where is the function defined for t G M by : 

= 1 r -,(*)- fU{p) < g{Y) < + y X r„) . 

|C(p)<9W<€(p)+2/xr„| V / 

For n sufficiently large, the function /i^(p) „(•) has Hermite rank Tp. Denote by cj^n the 
j-th Hermite coefficient of From Lemma El there exists a sequence {i^j)j>-p 

such that, as n ^ +00 

Since, for all n > 1 E(/i„(y)^) = J2j>Tpi'^j,n)'^ / J^- < +00, it is clear that the sequence 
(^j)j>rp is such that X]j>Tp('^i)^/i' < +00. By denoting A a positive constant, we get, 
as n — > +CXD 

E(W^„ - Zn? = ^ -^^^ E (/.^(p),„ {Y{^,)) /.^(p),„ (y(^2))) 

—2-1 " 

= ^ 7777-^^2 E E Ci„„C,„„E (y(n)) (Ffe))) 



n 



-2 T _ c2 



^ ^ ? E ^ E = ^ f ^ E P(^A = '^(^n), 

j>Tp \i\<n \ \i\<n j 

from Remark [TJ Therefore, Wn — Zn converges to in probability, as n — > +00. Thus, 
for all e > 0, we have, as n — > +(X), 

F(K <y,W^„ >y + e) =P(^n <y„,vr„ >y + e) ^0. 
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Following the sketch of this proof, we also have P (Ki > y + £, W„ < y) — 0, ensuring 
condition (b). Therefore, W„ — Z„ converges to in probability, as n — +oo. Thus, for 
all e > 0, we have, as n — ^ +oo, 

Following the sketch of this proof, we also have P (V^ > y + £, W„ < y) — 0, ensuring 
condition (b). 

Acknowledgement The author would like to thank the referee for his coments im- 
proving the statement of Corollary 2. 
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